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Abstract 

The theory of Kolmogorov-Johnson-Mehl-Avrami (KJMA) for phase 
transition kinetics is subjected to severe limitations concerning the 
functional form of the growth law. This paper is devoted to side step 
this drawback through the use of correlation function approach. More- 
over, we put forward an easy-to-handle formula, written in terms of 
the experimentally accessible actual extended volume fraction, which 
is found to match several types of growths. Computer simulations 
have been done for corroborating the theoretical approach. 
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1 Introduction 



The Kolmogorov- Johnson-Mehl-Avrami (KJMA) model [TJ |2J [3] finds appli- 
cation in a vast ambit of scientific fields which ranges from Thin Film Growth 
to Materials Science [U El El El [TOl [11] to Biology and Pharmacology 
[T2l IT3] , let alone the Applied Probability Theory [TJ1 [15]. In the majority 
of these studies the authors made use of a simplified version of the KJMA 
formula: the stretched exponential X(t) = 1 — exp(— at n ), where X is the 
fraction of the transformed phase, a and n (the latter known as Avrami's ex- 
ponent) being constants. The model, in principle, is simple because it rests 
on a Poissonian stochastic process of points in space, to which a growth law is 
attached. In fact, owing to the Poissonian process, the nucleation takes place 
everywhere in the space i.e. also in the already transformed phase. This par- 
tially fictitious nucleation rate (I(t)), for we are dealing with a Poissonian 
process, is linked to the actual (real) nucleation rate (I a (t)) according to: 
I a (t) = I(t)[l — X(t)], where X(t) is the transformed fraction. The growth 
law transforms each point in a nucleus of radius R(t), t stands for time. The 
pair, "points' generation" and "growth law", is a key quantity of the theory. 
It happens that the KJMA model fails for time-dependent points genera- 
tion rate (i.e. nucleation rate) associated with diffusional-type growth laws 
[16] , [17] . In particular, let us define two classes of growth laws: i) d 2 Rj dt 2 > 
and ii) d 2 R/dt 2 < 0. The KJMA model is suitable for describing the first 
class of growths and for this reason we named it KJMA-compliant as opposite 
to the second to which we attach the adjective KJMA-non-compliant. The 
reason for that is due to the particular stochastic process taken into account. 
As a matter of fact, the Poissonian process requires that points can be gen- 
erated everywhere throughout the space independently of whether the space 
is, because of growth, already transformed or not. Points generated in the 
already transformed space are named phantoms after Avrami. It goes with- 
out saying that, in the case of KJMA-compliant growths, phantoms do not 
contribute to the true transformed fraction, i.e. they are just virtual points 
whose only role is to simplify the mathematics [18J. On the other hand, in the 
KJMA-non-compliant growths phantoms may contribute to the phase tran- 
sition through the non-physical "overgrowth" events [3J. Incidentally, it is 
worth noticing that the KJMA-non-compliant growths and KJMA-compliant 
growths are indistinguishable if associated with simultaneous nucleation [19] . 

According to what has been said, one can summarize saying that the 
concept of phantom implies the existence of the two classes of growths. 

One has to bring in mind that the above stretched exponential expression 
( X(t) = 1 — exp(— at n )) is the exact solution of the kinetics only in the case 
I(t) = constant and I(t) ~ S(t), I being the nucleation rate and 8(t) Dirac's 
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delta function, provided the growth is according to a power law. In general, 
the term at n is a simple way to approximate the convolution product be- 
tween the nucleation rate (phantom included) and the nucleus volume. This 
convolution is the "extended" transformed fraction and takes into account 
the contribution of phantoms, X e . In view of the large use of the KJMA 
theory for dealing with experimental data, it should be desirable to make 
use of an extended transformed fraction deprived of phantom contribution, 
X e . A pictorial view of the geometrical meaning of X e , X e and X is reported 
in fig.l. 




Fig.l 



Figure 1: Pictorial view of a random ensemble of actual (dots 1-3) and phan- 
tom (dots 4-5) nuclei. The definitions of both phantom included (X e ) and 
actual (X e ) extended fractions are also reported. The (extended) volume of 
the i-th nucleus is denoted as Xi and the transformed fraction as X. 

The aim of this contribution is twofold: i) to model the phase transfor- 
mation kinetics in terms of actual quantities such as the nucleation rate; ii) 
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to provide an expression for the transformed fraction as a function of X e . 

2 Theory 

In this section we discuss the stochastic theory that has to be employed in 
order to get rid of phantoms in modeling the kinetics of phase transformations 
ruled by nucleation and growth. To this end, let us define the phantom 
included nucleation rate I(t), and the "actual" nucleation rate, I a (t), namely 
the rate of birth of the "real" nuclei; I a {t) is the quantity experimentally 
accessible. It goes without saying that a mathematical formulation of the 
phase transition kinetics which employs the "actual" nucleation rate holds 
true for both KJMA-compliant and KJMA-non-compliant growth laws. This 
automatically overcomes the limit of the KJMA approach. However, actual 
nuclei imply a severe complication of the stochastic-nature of the process 
under study, for we shift from a Poissonian to a non-Poissonian process: in 
fact actual nuclei are spatially correlated. Let us address this point in more 
detail by denoting with R(t, t') the radius of a nucleus, at running time t, 
which starts growing at time t' < t. To be an actual nucleus it has to lie 
at a distance r > R(t',t") from any other "older" nucleus with t" < t' . In 
other words, in the spirit of the statistical mechanics of hard spheres, this 
condition is formalized (at the lowest order) through the pair distribution 
function for the "pair of nuclei (t',t")" at relative distance r 

f 2 (r,t',t")nH(r-R(t>,t")), (1) 

where H(x) is the Heavyside function. Throughout the paper we employ 
the notation by Van Kampen [20] according to which n— dots distribution 
and correlation functions are denoted as /„ and g n , respectively. 

In previous papers, we have presented a theory for describing phase tran- 
sitions in the case of spatially correlated nuclei and for time dependent nu- 
cleation rate [2T1 122] . The untransformed fraction can be expressed in terms 
of either the distribution functions, (f n - functions), or the correlation func- 
tions, (^-functions), where the / n 's and g n 's are linked by cluster expansion 
|20j . Given a generic point of space we have computed the probability that 
this point is not covered (transformed) by any nucleus up to time t. This 
probability is the fraction of untransformed phase i.e. Q(t) = 1 — X(t). By 
denoting with A(t,t') the volume of a nucleus, which starts growing at time 
t' < t, at running time t, and in the case of symmetric f n and g n functions, 
the uncovered fraction is given by 
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Q(t) = 1-/ I(h)dtx I fi{vi)dv 1 + (2) 

JO J A(t,t{) 

+ I /(ti)tfti / I(t 2 )dt 2 f tin / dr 2 / 2 ( ri ,r 2 )-.... 

JO JO J A{t,t x ) J A(t,t 2 ) 

= 1 + 

+ V f- / /(ti)dti.. / /(tm)d* m X 

m=l m - 70 

/ dr x / dr 2 .... / / m (ri, ..,r m )cfr m 

JA(t,ti) JA(t,t 2 ) JA(t,t m ) 



or, 



g(t) = exp[^^ — >— I Htjdh... I(t m )dt m 

m=l m ' ^° ^° 

/ dr x / dr 2 .... / g- m (r 1; .., r m )(ir m . (3) 

JA(t,h) JA(t,t 2 ) J A(t,t m ) J 

It is worth pointing out that the nucleation rates entering these equations 
are in fact subjected to the condition imposed by the correlation among 
nuclei. This quantity may or may not imply phantoms depending on the 
specific form of the f n functions. For this reason we introduce the new 
symbol J. In particular, for the hard core correlation (eqnjl) J coincides 
with the actual nucleation rate, I = I a , which leads to the solution of the 
phase transition kinetics in terms of the actual nucleation rate. 

Since eqn,2]3 are the exact solutions of the stochastic process linked to 



the phase transition, they also coincide with the KJMA formula provided the 
above mentioned preconditions are met i.e. random nucleation and KJMA- 
compliant growth. 

As far as the hard-core correlation (eqnjl]) and the kinetics eqns, 2p 
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concerned, we note that the number of nuclei of size R(t,ti) is I a {ti)dt\ = 
0(dti). As a consequence and in the framework of the statistical mechanics 
of hard sphere fluid, we are dealing with an extremely dilute solution of pairs 
of components "ii,t 2 ". Accordingly, the number density of £ 2 spheres being 
of the order of 0(c?t 2 ), higher order terms in the cluster expansion of the / 2 
function can be neglected, thus f 2 {r,t',t") = H(r - R(if,t")) in eqn(2) 

2.1 KJMA-compliant growths 

In this section we show that eqn|2] is compatible with the KJMA kinetics 
only in the case of KJMA-compliant function. On the other hand, such 
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a comparison will give a deeper insight into the reasons why the KJMA 
kinetics does not work in the case of "KJMA- non-compliant" functions. In 
the following, we discuss the linear growth law for 2D case. Also, to simplify 
the complexity of the computation the "actual" nucleation rate, I a , is taken 
as constant; as a consequence the phantom included nucleation rate reads 
[31 E9] /(£) = / a /(l - X{t)) = I a /Q{t) = I a F{t) and the KJMA kinetics 
becomes 

F(t) = exp\ [ I a F(t')7rR 2 (t,t')dt'], (4) 
L Jo J 

where R(t, t') = v(t — t') and v is a constant. 

We consider the series expansion of F(t) around t = 0. One gets, 

_ F (n + 1) = {m){ n) = ^ ( n )F (n) nM) (g) 

fc=0 

where fi(t) = 27r/ a /J F(t')R(t, t')d t R(t, t')dt' = 2itl a v /* F(t')R(t,t')dt' . 

Moreover, since fi( m+2 ) = reF^, being k = 2irv 2 I a and = = 
fiW = 0, from eqn|5| it is found that only the terms F^ 3n - ) (0) are different 
from zero, i.e. F(t) = ^^° =0 j^F (3n \0)t 3n = ^2Z=o c mT ■ In particular, the 

first coefficients are: F( 3 )(0) = k,F^(0) = 11k 2 ,F^(0) = 375k 3 , F^ 12 \0) = 
234147k 4 . 

Next we derive the series expansion of the untransformed fraction Q(t) = 
}~2n=o b nt n , by exploiting the condition 1 = Q(t)F(t) = Y,n=o M n Em=o °mt m - 
Even in this case the b n coefficients are different from zero for n = 3k (with 
integer k). The first four coefficients are 



b = ~ = 1 

63 = ~° 3 2 (6) 

&6 = —Cq + C 3 

bg = —Cg + 2C3C6 — C3. 

By using the expression of c n 's the series of the untransformed fraction 
up to t 9 is given as 

Q(t) = 1 - iKt 3 + ^ 2 t 6 - Sf K 3 t 9 + 0(t 12 ). (7) 



Since X e = ^y/tt 3 = |7rJ a t> 2 t 3 = J ' I a %R 2 (t,t')dt' the series eqn 
rewritten as 

Q 6Q 
Q(t) = 1 - X e + —X 2 - — A 3 + 0(A e 4 ), 
vw 20 e 560 e V eJ 



can be 
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the two last coefficients being 0.45 and 0.123. We emphasize that the 
coefficients of this series only depend upon growth law, for constant I a (see 
also the last section). 

The next step is to show that the untransformed fraction, given by the 
series eqn|2]is equal to eqn|7]or eqnjSJ We have carried out the first two terms, 
exactly, while, owing to the tremendous computational complexity, for the 
third term an approximation has been employed. It is worth reminding that 
the distribution functions are /i(ri) = 1, f 2 (ri 2 , ti, h) = H(r 12 — R(ti,t 2 )) 
and /3(ri 2 ,ri3,r 2 3,ii,t2,*3) = H(r 12 - R(t 1 ,t 2 ))H(r 13 - t 3 ))H(r 23 - 
R(t2,t 3 )), where ry is the relative distance. In fact, since l a dti = 0(dti) 
the system of dots is dilute and f 3 (l,2,3) = / 2 (1, 2)f 2 (l, 3)^2(3, 2) i.e. the 
superposition principle holds true [23J . Since the system is homogeneous the 
f n functions depend on ry = |r^-| = — Eqn|2] becomes 

Q(t) = l-I a J*dtif A(tM) dr 1 + 

+ J a Jo dt ± Jo' df 2 / A (t,ti) dVl L{tM) rfl "2#( r 12 - R(h,t 2 )) + 

(9) 

~ J a Jo dt ± Jo 1 dt2 dt 3 lA(t,ti) dTl /A(t,t 2 ) ^ X 

f A{tts) dr 3 H{r 12 - R{t u t 2 ))H{r 23 - R(t 2 , t 3 ))H(r 13 - R(t u t 3 )) + 

where the integration domain A(t,tj) is the circle of radius R(t,ti) = 
v{t — ti). The fi containing term is the extended surface fraction X e (t) = 
I a f Q dti\A(t, ti)\ and coincides with the second term of the expansion eqn 8 
Let us focus our attention on the integrals in the spatial domain - for the 
sake of clarity shown in fig. 2a together with the circle of correlation R(t\,t 2 )- 
of the f 2 containing term. By employing relative coordinates the integrals 
read 



/ dn / dr 12 H{r 12 -R{t 1 ,t 2 )) = / dx x A{r x , t, h, t 2 ), (10) 

^A(t,ti) J A(t,t 2 ) ^A(t,<i) 

where A(ri,t,ti,t 2 ) = A(ri, R(t, t 2 ), R(ti, t 2 )) is the area "spanned" by 
the second nucleus when the first one is located at ri (t 2 < ti < t). It 
is at this point of the computation that the growth law comes into play; 
indeed in the case of KJMA-compliant growth laws (here linear growth) the 
correlation circle R(t\,t 2 ) is entirely within the integration domain R(t,t 2 ) 
(fig. 2a). Consequently, 
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Figure 2: The integration domains and the correlation circles are depicted 
for the integrals over the f 2 function. The cases of KJMA-compliant and 
KJMA-non-compliant growths are reported in panels a and b, respectively. 
In the drawing R u = R(ti,t 2 ), Ri = R(t,ti) (i = 1,2) and t 2 < t\. In the 
case of KJMA-compliant growth R\ + R\ 2 < R 2 and the hard core circle is 
within the integration domain, R 2 (a). In the case of KJMA-non-compliant 
growths R\ + R\ 2 > R 2 and the hard core disk overcomes the integration 
domain of the "older nucleus" (b). 
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A(t,t u t 2 ) = ir[R 2 (t,t 2 )-R 2 (t 1 ,t 2 )] (11) 

is independent of iv On the other hand, in the case of KJMA-non- 
compliant growth laws the correlation circle overcomes the integration do- 
main of the second nucleus and the relationship above does not hold true, 
anymore (see fig. 2b). In general, the KJMA-compliant functions satisfy 
the condition R(t,ti) + R(ti,t 2 ) < R(t,t 2 ) i.e. for a power growth law 
(t — ti) n + (t\ — t 2 ) n < (t — t 2 ) T \ which is verified only for n > 1. In 
fact, by setting r = t\ — t 2 and rj = t - z ^ 2L > 1, the inequality above reads 
(77 — l) n < 7] n — 1, which is satisfied for n > 1. On the other hand, for 
n = l/k (with integer k > 1) the inequality is [(77 — l) 1 '* + l] k < 77 namely, 

fc-i 

£(J)fo-l)*<0, (12) 
which is never satisfied (77 > 1). 

For the KJMA-compliant growth the contribution of the f 2 containing 
term becomes 



wi 2 a Jj dh f* 1 dt 2 J A(Mi) dn [/? 2 (t, t 2 ) - i?(t X , t 2 

= T a J* Jo Jo 1 dt * R2 ^ *0 i R2 ^ *a) - ^1^2)] (13) 



*2 
2 



It is possible to show that for linear growth the last term of eqn 13 is equal 
to j^(irl a v 2 t 3 ) 2 , consequently we get ^f- — = §qX 2 , which coincides with 
the term of the same order in the KJMA series eqn|8j 

Let's now briefly consider the contribution of the f 3 containing term in the 
general expression eqn|2] Fig. 3 shows the integration domains Ri = R(t, tj) 
and the correlation circles R(U,tj) for the three nuclei born at time £j (i = 
1,2,3). The configuration integral over f 3 in eqn|9] becomes 



J a lo dti /J 1 dt 2 J ' 2 dt 3 f A(tM) drt / A(t>ta) dr 2 H(r l2 - R(t x , t 2 ))x 
J A(t t3) dr 3 H(r 13 - R(t x , t 3 ))H(r 23 - R{t 2 , t 3 )) 

(14) 

= II ft dt, £ dt 2 dt 3 x 

Ia^i) dvi lA{t,t 2 ) dr 2H(r 12 - R(t u t 2 ))A(r 12 , R(t u t 3 ), R(t 2 , t 3 )), 
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0<t 3 <t 2 <t 1 <t 























Fig.3 



Figure 3: Integration domains and correlation hard disks for the integral 
over the function. In the drawing R^j = R(ti,tj) and Ri = R(t,ti). In the 
case of KJMA-compliant growth all the circles are entirely within the R(t, t^) 
circle. 
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where £3 < £ 2 < £1 < £ is assumed. In this equation, A(ri 2 , R(t\, £3), R{t 2 ,t 3 )) 
is the area spanned by the third nucleus when the first and the second are 
located at r x and r 2 , respectively. Because of the possible overlap between 
the correlation circles R(ti, £3) and R(t 2 , £3), and for R(t±, £3) is encompassed 
within R(t, £3) this area is a function of r\ 2 as 

A(n 2 ) = u{r 12 )H{R{t u t 3 ) + i2(t 2 , £ 3 ) - r 12 ) + 

(15) 

+7T [i? 2 (£, £ 3 ) - R 2 (h, £ 3 ) - i? 2 (£ 2 , £3)] #(r i2 - i?(£ 2 , £3) - R(h, £3)), 

where w(r 12 ) = 7r(i? 2 (£, t 3 )-R 2 (t u t 3 )-R 2 (t 2 , t 3 ))+w(r 12 , R(h,t 3 ), R(t 2 ,t 3 )), 
with zu(x, pi, p 2 ) being the overlap area of two circles of radius pi and p 2 at 
relative distance x 



w(x,p u p 2 ) = -\^/Ax 2 p\ - \p\ - x 2 - p 2 } 2 + 

(16) 

+pf arccos 1 2xpi + p 2 arccos 2 2rrpa ■ 
It turns out that the computation of the third order term of the series, 



eqn,14, is a formidable task indeed. We do not attempt to perform the exact 
estimate of this term which, however, must coincide with the same order 
term of the KJMA series. On the other hand, an approximate evaluation of 
this term by using an oversimplified form of the A(ri 2 ) area, is possible by 
formally rewriting this area as A(r 12 ) = 7r[i? 2 (£, £ 3 ) — R 2 (t%, £3) — (3R 2 (t 2 , £3)], 
where (3 e (0, 1) is given by = [1 - ^g^H{R{t u t 3 ) + R{h,t 3 ) - r 12 )]. 
In the case of complete-overlap (/3 = 0) we get, 



tt 3 /3 f*dti C dt 2 f* 2 dt 3 R 2 (tM)[R 2 (tM) - i? 2 (£i,£ 2 )p 2 (£,£ 3 ) - R 2 (ti,t 3 )], 

(17) 

that for the linear growth (R(t,f) = v(t - £')) gives ^I^n 3 vH 9 = ||X e 3 to 
be compared with the exact value ^ X 3 which brings an uncertainty of 4.3% 



2.2 KJMA- non-compliant growths 

Let us now consider the parabolic growth R(t, £') = v \/t — £'. In this case 
the series expansion of the function F(£) = 1/Q(t), given by eqnjij can be 
performed by employing the same computation pathway discussed above, 
where now = kF^ 1 ' and k = nl a v 2 . In this case F^ 2n \0) implying 

F(t) = 1 + f£ 2 + \n 2 t A + f kH & and 
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Q(t) = l-p 2 + £t 



180 



t 6 + 0(t s 



-,xi 



!IO-'e + 



the two last coefficients being 0.33 and 0.044. It is worth pointing out 
that in such an evaluation the transformed fraction, X, is comprehensive of 
the contribution of phantoms. In fact, we recall that eqn|4] is the KJMA 
solution with the phantom included nucleation rate, I a /(1 — X). 

The fx containing term of eqn|2] gives the extended volume fraction |t 2 . 
As far as the third term is concerned (f 2 contribution), it is possible to 
show that also for n = 1/2 the integral eqn,13 coincides with the third term 

does not 



Y^t 4 . However, it is important to stress that eqn 
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of eqn 

coincide, in this case, with the integral over the f 2 function of the exact 
solution eqn|9j since in the latter equation enter the actual nuclei, only. 
From the mathematical point of view, in the case of parabolic growth, the 
correlation circle is not contained within the integration domain as depicted 
in Fig. 2b. In other words, for KJMA-non-compliant growth laws the area A 
is a function of r\ and the term of order I 2 in eqn|9] does not coincide with 
(parabolic growth). In particular, under these circumstances 



-f- of eqn 
we get 



A(r u t, h, t 2 ) = u{n)H{n + R(h,t 2 ) - R(t, t 2 ))+ 



(19) 

+7T [R 2 (t, t 2 ) - R 2 {h,t 2 )] H(R(t, t 2 ) - n - R{h, t 2 )), 
where oj{r\) = nR 2 (t, t 2 ) —w(rx, R(t, t 2 ), R(t\, t 2 )) with zu(r\) the overlap 



area of two circles of radius R(t,t 2 ) and R(ti,t 2 ) at distance r\ (eqn. 16). 



3 Numerical Simulations 

The ultimate aim of this section is to propose a simple formula for describing 
the kinetics on the basis of the "actual" extended transformed fraction, X e . 
On the ground of eqn|3] the transformed fraction can be rewritten in the 
general form 

X(t) = 1 - exp[-X e (t) 7 (X e (t))], (20) 

where j(X e (t)) embodies the contributions of correlations among nuclei 
|21j . It is worth noticing that for KJMA-complaint growths (random nucle- 
ation) eqn {3] actually reproduces the KJMA formula. In fact by identifying / 
with the phantom included nucleation rate (1 = 1) one gets g m> \ = leading 



12 



to the formula Q = expf— X e ). On the other hand, working with the actual 
nucleation rate, in eqn|3j I = I a , g m ^ and the series has infinite terms. 
7(A" e (t)) can therefore be expanded as a power series of the extended actual 
volume fraction, X e . Moreover, by exploiting the homogeneity properties of 
the f n functions (see below), it is possible to attach a physical meaning to the 
power series coefficients, in terms of nucleation rate and growth law. Also, 
for constant I a the coefficients of this series only depend upon growth law. 
To the aim of achieving a suitable compromise between handiness and plia- 
bility, we retain the liner approximation of ~/(X e (t)) obtaining the following 
kinetics 



X = 1 -exp[-(aX e + &X e 2 )], (21) 

with a and b constants. For the sake of completeness we point out that, 
according to its physical meaning, the parameter a should be unitary. Nev- 
ertheless, the substitution of the infinite expansion with only two terms au- 
thorizes the introduction of the new parameter a. In any case, the a values 
are found to be nearly one (see fig. 7 below). 

In order to study the transition kinetics in terms of actual nuclei and to 
test eqn,21, we worked out 2D computer simulations for several growth laws 
at constant nucleation rate, I a . 

As typical for this kind of study [5], [21], the simulation is performed on 
a lattice (square in our case) where, in order to mimic the continuum case, 
the lattice space is much lower than the mean size of nuclei. 

In particular, the transformation takes place on a square lattice whose 
dimension is 1000 x 1000 with a nucleation rate of I a = 3. It is worth 
reminding that, since the nucleation is Poissonian, it occurs on the entire 
lattice independently of whether the space is already transformed or not. 
The computer simulation can be run taking into account the presence of 
phantoms or not. In the former case the outputs have been labeled as "w" 
while the latter as "wo". As far as the growth laws are concerned, we limited 
ourself to the power laws, R(t) ~ t n for n = 1/4, 1/3, 1/2, 1, 3/2 and 2. 

The results of the simulations are displayed in figs.4a-c for the KJMA- 
non-compliant growths. In particular, the fractional surface coverage,A, as 
a function of the actual extended fraction X e , with and without the con- 
tribution of phantoms, are reported (curves labeled with "w" and "wo", 
respectively). 

The contribution of phantom overgrowth to the transformation kinetics 
is highlighted in fig.5 and shows that this effect brings an uncertainty on X 
which ranges from 2% to 5% on going from n = 1/2 to n = 1/4. In the case 
of parabolic growth this figure is lower than 2%. 
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Figure 4: Computer simulatiosn of phase transformations ruled by KJMA- 
non-compliant growths. The surface fraction, X, is shown as a function of 
the extended fraction X e , for the power law (R « t n ) where n = 1/2, 1/3 
and 1/4 in figures a, b and c, respectively. The kinetics with (w) and without 
(wo) the inclusion of phantoms are displayed. 
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Figure 5: Contribution of phantom overgrowth to the kinetics reported in 
fig. 4 (displayed is the difference between curves (w) and (wo)). The area 
beneath the curves normalized to the area of the kinetics of fig.4 are, respec- 
tively, 0.053, 0.038 and 0.02 for n = 1/4, n = 1/3 and n = 1/2. 



These results are in qualitative agreement with previous studies on phan- 
tom overgrowth, although performed for a different nucleation laws [2H HE] . 
As discussed in more details below, the results displayed in figs. 4, 5 are uni- 
versal, i.e. they only depend on power exponent, n, and nucleation law (in 
the present case I a = constant). Accordingly, the lower n the more im- 
portant is phantom overgrowth. In fact, let us consider a phantom, which 
starts growing at time i, located at r p from the center of an actual nucleus 
which starts growing at t = (fig.6a). For KJMA-non-compliant growth, 
R(t — t') = v(t — t')k (with integer k > 1), the phantom overtakes the actual 

nucleus at time t Q , that is the solution of the equation r p + v(t a —t)k— vt£ , 
namely 



£ = fa- 1)*, (22) 



where £ = -^r < 1 and r] — y > 1. The graphical solution of eqn 

Vt Ti 
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is 
k. 



depicted in fig. 6b and indicates that t a (and therefore rj) decreases with" 
This is in agreement with the results of fig. 5 which shows that the overgrowth 
phenomenon is more important at greater k. 



As far as the guess function eqn, 21 is concerned, it matches the simulation 



curves with a very high degree of correlation. For instance the output of the 
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Figure 6: Sketch of the overgrowth process in the space domain. r p denotes 
the location of the phantom which start growing at time t when the size of 
the actual nucleus is R{t). The phantom overtakes the actual nucleus at time 
t Q when the size of the actual nucleus is R(t ). 
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fit to the n = 2 curve, gives a = 0.9750 ± 0.0007, b = 0.088 ± 0.001 and a 
squared correlation coefficients practically 1. For the sake of completeness 
the a and b fitting parameters are shown in fig. 7, where a is found to be nearly 
one. This is in agreement with the theoretical value predicted by eqnJHJ 



0.1 



0.01 




0.4 0.6 0.8 1 

n 



Figure 7: Behavior of the fitting parameters a and b of eqn,21 as a function 
of growth exponent n, for KJMA-non-compliant growths. 



The behavior of the transformed fraction for KJMA-compliant growths 
are reported in fig.8 for n = 1, n = 3/2 and n = 2. 

These kinetics are very close to each other and differs, markedly, from 
that at n — 1/2 also reported in the same figure. In the inset, the kinetics 



for n = 1/2 and n = 1 are compared with the KJMA series expansion eqn.18 



and eqn|8j respectively. The fact that the curves for n = 1, 3/2 and 2 collapse 
on the same curve, can be rationalized computing the coefficients of the series 
expansion of Q(X e ) for integer n. In particular, by employing the method 
discussed in the previous section the two last coefficients of the series (e.g. 
eqn|8]) are 0.4960, 0.1633 and 0.5, 0.1673 for n = 2 and 3, respectively. We 
also performed computer simulations of phase transitions for non-constant 
actual nucleation rate. The output of this computation is displayed in fig. 9 
where the behavior of the nucleation rate is shown in the inset as function 
of X e . In particular, the actual nucleation rate is given by the function 
I a {t) w t 2 exp(-at 3 ). 

Also in this case the function eqn,21 has been found to match the kinetics 
with high degree of correlation where, again, the independent variable is the 
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Figure 8: Kinetics of the actual surface fraction, as a function of actual 
extended surface fraction, for several values of n. In the graph the kinetics 
for n = 1/2 is compared with the KJMA-compliant growth with n = 1, 
n = 2 and n = 3/2 (from the top, respectively). In the inset the kinetics for 
n — 1 and n = 1/2 are displayed together with the truncated KJMA series 



expansions eqn|3 and eqn 18 , respectively 
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Figure 9: Fig. 9 Kinetics of the actual surface fraction, as a function of the 
actual extended surface, for non-constant actual nucleation rate (left scale). 



The best fit of the function eqn 21 to the X(X e ) kinetics has been shown as 



dashed line. The correlation coefficient of the fit is in fact 1 for the parameters 
a = 1.014 and b = 0.0382. The actual nucleation rate, as a function of X e , 
is also reported on the right scale in nuclei x 10 6 / site units. 
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actual extended surface fraction. 

Let us address in more detail the question of the dependence of volume 
fraction on extended volume fraction. To this end, we discuss the second 
order term of the exact solution eqn|2j namely 

f I a (t')dt' f Ia(t")dt" [ dr x A{r X iR(t'it"),R{t,t")) } (23) 

Jo Jo J &(t,t') 

where eqn, 10 has been employed. We point out that in eqn|23|A(ri, R(t', t"), R(t, t")) 
is a second order homogeneous function of r±, R(t', t") and R(t, t") variables. 
Accordingly, for the growth law R(t,t') = v(t — t') n , using the re-scaled 
variables, r[ = ri/vt n , r' = i! jt and r" = t"/t the integral becomes 

K 2 vH 4n + 2 f 1 Ia{r , )dr , f T Ia{j » )dT » [ dr ' lA (r' 1} r'y'). (24) 



JO JO JA(t') 

Eqn, 24 takes the form C2X 2 where 

n , r T f m Jo dr'I a {r') I a (r")dr" f , A(r[S ,r")dr[ 

C 2 (n, [Ia(T)\) = $ i a (T>)(i-T>)*>dT>)2 depends on n and ac- 

tual nucleation rate. It is apparent that in the case discussed in the previous 
section I a (t) = I a , ^(n), as well as higher order coefficients, is a function 
of n, only. In this case the transformed fraction is expected of the form 
l-X = Y, k C k {n)Xt 

On the other hand, in the case of a constant "phantom-included" nucle- 
ation rate, I a (t) = Iq(1 — X(t)) = IoQ(t), and eqn J2] becomes an integral 
equation for the Q(t) unknown. With reference to the second order term, 



in this case eqn 24 takes the general form C'(n)X^ which now implies the 
series 1 — X = J2k ^'k( n )Xe (note that this is a series expansion in terms of 
phantom-included extended surface). In the specific case of KJMA-compliant 
growths, however, these series reduces to the exponential series with constant 
coefficients v , / . It is instructive to estimate the first two coefficients in the 
case of linear growth. For constant phantom included nucleation rate the 
untransformed fraction satisfies the integral equation, 



Q(t) = l-I J*Q(t>)\A(t,t')\dt>+ 

(25) 

+P Jl dt> J*J dt"Q(t')Q(t") J A(t tl) drUfat, t', t") + O(/ 3 ). 



The first order term of this equation, namely of the order of Iq 
Q(t) = 1 — I nv 2 t 3 /3 = 1 — X e . By substituting Q rs 1 — X e in eqn 
get 



, gives 



25, we 
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Q(t) = 1-X e + I f*X e {t')\A{t,t')\dt'+ 

(26) 



+/ 2 / * dt> ft dt" J dv.Ain, t, t', t") + o(ff 



Using dimensionless variables r[ = ri/vt, r' = t'/t and r" = t"/t eqn.26 
eventually becomes 

Q(t) = 1-X e + 3X 2 e £ r' 3 (l - r') 2 dr'+ 

(27) 



+9X e 2 £ dr' / T dr"(l - r'fA{r', r") + 



where A(r', r") is given throug h eq njllj Notably, the last term in eqn,27 
has been already estimated in eqn.13 and is equal to ^X 2 . The coefficient 



of X 2 is eventually computed as ^ + ^ = |, that is the expected result. 

It is worth noting that the present approach can also be applied to differ- 
ent convex shapes other than circles and spheres, provided the orientation of 
nuclei is the same (with a possible exception for triangle) . This aspect has 
been discussed in details in refs.[25].[26]. 

We conclude this section by quoting the recent results of ref. [17J. In 
this noteworthy contribution the author faced the problem of describing the 
kinetics in terms of the actual nucleation rate. An ingenious application of 
the so called "Differential critical region" approach makes it possible to find 
the Q(t) kinetics by solving an appropriate integral equation [T7j. On the 
other hand, the different method employed in the present work, based on the 
use of correlation function, pertains to the same class of stochastic approaches 
on which "Kolmogorov's method" is rooted. It could be enlightening, for the 
present topic, to demonstrate that the two approaches are in fact equivalent. 



4 Conclusions 

We have shown that, employing the correlation function approach, the con- 
straints on growth laws underlying the KJMA theory can be eliminated. In 
other words, the present modeling is not constrained to any form of the 
growth law. The actual extended volume fraction is shown to be the natural 
variable of the kinetics, which implies universal curves. Besides, we proposed 
a formula to fit experimental data by using the measurable actual extended 
coverage. The displacement of the kinetics from the exponential law, i.e. the 
b parameter in eqn|2T| may give insights into the microscopic growth law of 
nuclei. 
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